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' There are three parameter degeneracies inherent in the three-neutrino 

■ analysis of long-baseline neutrino experiments. We develop a systematic 

method for determining whether or not a set of measurements in neutrino 
Q-i! oscillation appearance experiments with approximately monoenergetic beams 

Q_i! can completely resolve these ambiguities. We then use this method to identify 

experimental scenarios in which the parameter degeneracies may be efficiently 
resolved. Generally speaking, with two appearance measurements degenera- 
cies can occur over wide areas of the (6,613) parameter space; with three 
measurements they occur along lines in the parameter space and with four 
measurements they occur only at isolated points. If two detectors are placed 
at the same distance from the source but at different locations with respect 
to the main axis of the beam (a detector cluster) , each detector will measure 
neutrinos at different energies. Then one run with neutrinos and one run with 
antineutrinos will give the four independent measurements that in principle 
can resolve all of the parameter degeneracies if sin 2 2#i3 > 0.002. We also 
examine scenarios with detector clusters using only neutrino beams. Without 
detector clusters, the measurement of neutrinos and antineutrinos at a short 
distance and only neutrinos at a longer distance may also work. 



I. INTRODUCTION 



The recent solar neutrino data from the Sudbury Neutrino Observatory |T],0], which infers 
different neutrino fluxes from the charged- and neutral-current measurements, provides 
convincing evidence that electron neutrinos do in fact change flavor as they travel from 
the Sun to the Earth. The neutral-current measurement is also consistent with the 
solar neutrino flux predicted in the Standard Solar Model ||]. The preponderance of solar 
neutrino data increasingly prefers the Large Mixing Angle (LMA) solution to the solar 
neutrino puzzle, with bm\ x ~5x 10~ 5 eV 2 and amplitude close to 0.8 [§,fH|. This solution 
will be tested decisively by the KamLAND reactor neutrino experiment |J. 

The atmospheric neutrino deficit also gives a strong indication that neutrinos have mass 
and oscillate from one flavor to another - the most compelling interpretation is that u^s 
created in the atmosphere oscillate to v T with almost maximal amplitude and mass-squared 
difference Sm^ ~ 3 x 10~ 3 eV 2 [0. The K2K experiment || with a baseline of 250 km 
has preliminary results that are in agreement with this interpretation. Oscillations of 
to v e as an explanation of the atmospheric anomaly are ruled out by the CHOOZ || and 
Palo Verde |T(J reactor experiments, which place a bound on the — > v e oscillation am- 
plitude smaller than 0.1 at the 95% C.L. in the Sm^ region of interest. The MINOS ||11| , 
ICARUS [|T2[ and OPERA JTB| experiments are expected to come online in 2005 and study 
aspects of the oscillations at the atmospheric scale |14j]. The low energy beam at MINOS 



will allow a very accurate determination of the leading oscillation parameters. ICARUS 
and OPERA should provide concrete evidence that — > v T oscillations are responsible for 
the atmospheric neutrino deficit by identifying tau neutrino events. These long-baseline 
experiments, when combined with KamLAND, should provide further information on all of 
the parameters in the three-neutrino mixing matrix, except for the CP-violating phase (5) 
and possibly not the mixing angle associated with — > v e oscillations in atmospheric and 
long-baseline experiments (#13). It will take a new generation of experiments to provide 
accurate measurements of these parameters. 

Future measurements of #13 and 5 are not completely straightforward. The parameter #13 
cannot be determined from a single v measurement since the measured oscillation rate also 
depends on the value of 5, which is not known. Even with both a v oscillation measurement 
and a v oscillation measurement there are several parameter degeneracies that enter the 
determination of the three-neutrino mixing matrix. These include the ambiguities associated 
with (i) the parameters in the XJ& (=sin6 l i3 e~ tS ) element [ 15— |TqH , (ii) the sign of the mass- 



squared difference responsible for the oscillations of atmospheric neutrinos (^m 2 ^) p6|-p0 



and (iii) the primary mixing angle for the oscillation of atmospheric neutrinos (#23 ) [|17| , |T9 



Each of these ambiguities can mix CP violating (CPV) and CP conserving (CPC) solutions. 
To resolve the ambiguities requires multiple measurements; the usual prescription involves 
making measurements with different beam energies, baseline lengths, and/or beam particles 
(neutrinos versus antineutrinos). A typical scenario requires measurements of neutrinos and 
antineutrinos, each at two different energies [J17|| , which would entail four separate runs, each 
of which would take years to complete. 

In this paper we discuss alternate possibilities that capitalize on the fact that the en- 
ergy spectrum of a conventional neutrino beam varies with the angle from the main beam 
axis [pl| |23|1 . The off-axis components of a neutrino beam have very narrow spectra for a 
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given off-axis angle. Viewed this way, a single neutrino beam is actually a continuum of 
approximately monoenergetic beams; the energy can be tuned simply by varying the angle 
with respect to the main beam axis. This has the advantage that experiments at more than 
one neutrino energy may be run simultaneously using a single beam line, simply by placing 
detectors at slightly different angles from the main axis of the beam (a detector cluster). 
One could also locate one detector on the beam axis and another off-axis. Furthermore, 
the high energy tail of off-axis components of the beam is suppressed, reducing the high- 
energy contribution to backgrounds. One disadvantage is that the neutrino flux decreases as 
the off-axis distance increases, which may give a practical limit for the choices of neutrino 
energies. 

Detailed studies have already been made for off-axis beams with single detec- 
tors [ff9pl|,p5f ) but how to resolve the parameter ambiguities has not been thoroughly 



addressed. In this paper we present some possible scenarios for resolving parameter de- 
generacies using detector clusters with off-axis neutrino beams. In particular we examine 
how well parameter degeneracies may be resolved with one run of a neutrino beam and one 
run of an antineutrino beam. We also discuss other, more speculative, possibilities involving 
detector clusters, such as multiple detectors and only neutrino beams. We compare these 
results to the ability of more conventional experimental setups (without detector clusters) 
to resolve parameter degeneracies. 

The organization of our paper is as follows. In Sec. [Q| we present oscillation probability 
formulas based on a constant density approximation that works well for the situations that 
we consider. In Sec. |T| we determine in a general way the number and type of measurements 
that are needed to resolve parameter degeneracies. Specific detector scenarios are presented 
and discussed in Sec. [IVl. Concluding remarks are made in Sec. M. 



II. OSCILLATION PROBABILITIES IN MATTER 

We work in the three-neutrino scenario. For oscillation studies, the neutrino mixing 
matrix U can be specified by 3 mixing angles (#23 , 6*12, 6 Vi ) and a CP- violating phase 5. We 
adopt the parametrization 

/ C13C12 C13S12 s 13 e~ lS \ 

U = — C23S12 — Si3S 2 3Ci2e i<5 C23C12 — ■Sl3 s 23' S 12e" S C13S23 , (1) 

1 iS iS I 

\ S23S12 — Si3C23Ci2e — S23C12 — Si3C23Si2e C13C23 / 

where Cjk = cos^ and Sjk = smOj^. For Majorana neutrinos there are two additional 
phases, but they do not affect oscillations. In the most general U, the % are restricted to 
the first quadrant, < 9^ < tt/2, with S in the range < 5 < 2n. We assume that z/ 3 
is the neutrino eigenstate that is separated from the other two, and that the sign of Sm^ 
can be either positive or negative, corresponding to the case where h> 3 is either above or 
below, respectively, the other two mass eigenstates. The magnitude of 5rn^ 1 determines 
the oscillation length of atmospheric neutrinos, while the magnitude of Srn^i determines 
the oscillation length of solar neutrinos, and thus l^m^l <C IJmjj-J. If we accept the likely 
conclusion that the solar solution is LMA H^H, then Sm^i > and we can restrict 9 12 to the 
range [0, vr/4]. It is known from reactor neutrino data that #13 is small, with sin 2 2# 13 < 0.1 
at the 95% C.L. ||. Thus a set of parameters that unambiguously spans the space is Sm^ 
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(magnitude and sign), Sm^x, sin 2 2#i 2 , sin #23, and sin 2 26 l i3; only the #23 angle can be below 
or above tt/4. 

To determine oscillation probabilities we use the constant density approximation and 



expand in terms of the small parameters #13 and Sm^ [Pq ,P7||, which has been shown to 
reproduce well the exact oscillation probabilities for E v > 0.5 GeV, 6> 13 not too large, and 
L < 4000 km [[17]]. Up to second order in a and #13 , the oscillation probabilities for Sm^ > 
and Sm^ > are 

P(*V->i/ e ) = \x f + yge l{A+S) \ 2 = x 2 f + 2xyf g(cos5 cos A - sin 5 sin A) + y 2 g 2 , (2) 
Pfa^v*) = \xf + yge i(A ' s) \ 2 = x 2 f 2 + 2xyfg(cosScosA + sm5smA)+y 2 g 2 , (3) 



respectively, where 



x = sin 6 23 sin 2#i 3 , (4) 

y = a cos 623 sin Idyi , (5) 

/,/ = sin((l T i)A)/(l T i), (6) 

5 = sin(iA)/i, (7) 



and 



A = \Sm 2 31 \L/AE u = 1.27|5m 2 1 /eV 2 |(L/km)/(^/GeV) , (8) 
A= \A/5m 2 31 \, (9) 
a = l^m^/^m^l . (10) 

In Eq. (0), A/2E U is the amplitude for coherent forward charged- current z/ e scattering on 
electrons, with 

A = 2y/2G F N e E v = 1.52 x 10" 4 eV 2 Y e p (g/ cm 3 )E u (GeV) , (11) 

and iV e is the electron number density, which is the product of the electron fraction Y e (x) 
and matter density p(x). In the Earth's crust and mantle the average matter density is 
typically 3-5 g/cm 3 and Y e ~ 0.5. In all of our calculations we use the average N e along the 
neutrino path, assuming the Preliminary Reference Earth Model |28| . 

The coefficients / and / differ due to matter effects (A ^ 0). The values of /, /, and g 
(scaled by \jET v j L to account for the dependence of neutrino cross section with energy and 
the flux with distance) are shown versus A in Fig. |1] for L = 300 km (JHF to Kamioka), 
730 km (Fermilab to Soudan or CERN to Gran Sasso), 1290 km (Fermilab to Homestake), 
1770 km (Fermilab to Carlsbad), and 2900 km (Fermilab to SLAC); the scaling factor is 
chosen this way since the oscillation probabilities are quadratic functions of /, /, and g. 

To obtain the probabilities for 5m\ x < 0, the transformations A — > —A, y — > — y and 
A — > —A (implying / <-> — / and g — > —g) can be applied to the probabilities in Eqs. (0) 
and (]!|) to give 

P(vn ~^ v e) — x 2 f 2 — 2xyfg(cos5 cos A + sin S sin A) + y 2 g 2 , (12) 
P(V(i — ► & e ) = x 2 f 2 — 2xyfg(cosS cos A — sin 5 sin A) + y 2 g 2 . (13) 
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In practice, neutrino beams are not monoenergetic, even for a narrow band beam or 
an off-axis component of a beam. The number of appearance events is N(v^ — > u e ) = 
J P{Vfj, — > f e ) &crdE u , where $ is the energy-dependent neutrino flux, a the interaction 
cross section, and P the oscillation probability from Eq. (FJ). Then for Sm^ > 0, N has the 
approximate form 

Nfa — > f e ) = Aia; 2 + v4 2 a; cos 5 + A 3 x sin 5 + A A , (14) 

where 

A x = J f 2 $adE u , 

A 2 = J 2yfg cos A$(j dP,, , 

A 3 = ~J 2y fg sin A$adE u , 

A 4 = J y 2 g 2 §odE u . (15) 

There are similar expressions for N(pn — > u e ), and for Sm^ < 0. In all cases, whether one is 
describing the number of events at a single energy or total events integrated over an energy 
spectrum, the result can be written approximately as a linear combination of x 2 , xcosd, and 
xsin<5. In fact the probability for constant density can be written as A cos 5 + Bsin6 + C 
without any approximation |29j . It is this generic property that we exploit in this paper, so 



many of the qualitative aspects of our analysis should hold for neutrino beams that are not 
monoenergetic if only the total number of events is being used. 



III. DETERMINING THE OSCILLATION PARAMETERS 
A. Parameter degeneracies with two measurements 

It is expected that 8ra\ x and sin 2 2#23 will be well-measured in survival experiments; if 
the solar solution is LMA, as it now seems to be 0,f|||, then Sm^i and sin 2 2#i2 will be well- 
determined by KamLAND 0. Thus it is the parameters #13 and 8 that will be primarily 
measured in appearance experiments, modulo uncertainties in the other parameters. 

A single measurement is not sufficient to determine #13 since the oscillation probability 
depends on both #13 and 5 (see, e.g., Eqs. @ and @). In fact the inferred measurement 
of sin 2 2#i3 could be uncertain by an order of magnitude with only one measurement |T7 . 
The usual approach for determining # 13 and S is a measurement of both P and P at one 
L and E v . Then # 13 and 5 can be determined from Eqs. (0) and (Bf). Unfortunately, there 



are usually two solutions with distinct #13 and 5 (the well-known (5, #13) ambiguity |L5]-[n|). 
The existence of this ambiguity can be understood as follows. If #13 is held fixed, then as 5 is 
varied an ellipse is traced out in (P, P) space; the eccentricity of the ellipse is determined by 
the oscillation argument A P7| , ppfl . Ellipses for two different values of 9\ S can intersect, so 
that a single point in (P, P) space can be obtained by two different sets of parameters (5, 613) ■ 
For a given (P, P) there are also often two solutions for Sm^ < from Eqs. ([12]) and (|13D 



(the sgn(5m§ 1 ) ambiguity). Furthermore, since only sin 2#23 is known from atmospheric 
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neutrino experiments (or from measuring survival), if #23 7^ tt/4, each solution discussed 
above also has two possible values of #23 (each with a different #13 and S), one with #23 < 7r /4 
and one with 6*23 > 7r/4 (the (6*23,7r/2 — #23) ambiguity). Thus, in general, there can be as 
much as an eight-fold degeneracy of solutions from a measurement of P and P at one L 
and E v |T7|], which has already required two experimental runs (one with neutrinos and one 
with antineutrinos) . Further measurements at different oscillation arguments A are needed 
to eliminate all of the ambiguities, requiring additional runs. 

A judicious choice of the L and E v can help reduce the effect of the ambiguities. For 
example, the (5, 6*13) ambiguity can be reduced to a (5, n — 5) ambiguity by choosing L/E u = 
(2n — 1)(410 km/GeV)(3 x 10~ 3 eV 2 /{Sm 2 ^, where n is an integer (which implies A = 
in — |)7r, the approximate location of the appearance oscillation peak). The parameter 
#13 is removed from the ambiguity (and therefore in principle determined) since the ellipse 
in (P,P) space collapses to a line and the lines for different #13 no longer intersect |l6Hl8|. 
Furthermore, if L is sufficiently large, the large matter effects separate the ellipses for 8m\i > 



and 5m^i < 0, so that they do not intersect, removing the sgn^m^) ambiguity ||16| , |i7 



However, measuring (P, P) at a single L and E v will still leave a (5, n — 5) ambiguity, and 
may have a #23 ambiguity if #23 7^ tt/4 (with a possible associated CPV/CPC confusion), 
so that additional measurements are needed to completely determine the parameters and 
resolve all of the ambiguities. 

The arguments above can be generalized to show that any two measurements of neutrino 
and/or antineutrino appearance probabilities, even if they are not at the same L and E„, are 
subject to the same degeneracy problem. Each of the neutrino or antineutrino probabilities 
in an appearance experiment (e.g., Eqs. (0) and (|3|), or ( |P2"D and (|T3D) is a linear combination 
of the parameters x 2 , xcosS, and xsinfi. The coefficients of these parameters involve /, /, 
g, and A, which depend only on L, E u , Sm^, and y. It is not hard to show that for two 
such measurements of P and/or P, if #13 is fixed and 5 is varied, an ellipse is traced out in 
the two-dimensional probability space. Ellipses with somewhat different values of #13 will 
overlap, so that there is always a (5, #13) ambiguity for any two appearance measurements. 
The only exception is the special situation where the ellipse collapses to a straight line, in 
which case lines for different #13 no longer overlap, but there is still an ambiguity involving 
5 since the two halves of the collapsed ellipse then overlap (e.g., the (8, n — 5) ambiguity for 
(P, P) when A = tt/2). Similarly, there is a #23 ambiguity if #23 7^ 7r/4, and there may be 
a sgn(5m3 1 ) ambiguity if the matter effect is not large enough to separate the ellipses for 
8m\i > and Sm^ < 0. Thus two measurements will always have an ambiguity involving 5 
(which may or may not also involve #13, depending on the value of the oscillation argument 
A), and may also have ambiguities involving sgn^m^) and #23. Ambiguities involving 
sgn(5m| 1 ) and #23 will have a CPV/CPC confusion; the ambiguity that involves 5 alone 
will have a CPV/CPC confusion unless A = mr/2. 

We emphasize that the degeneracies discussed here are exact, i.e., there are different sets 
of parameters that give identical predictions. Thus these ambiguities are present even in 
the limit of no experimental uncertainties. When such uncertainties are taken into account, 
nearly degenerate solutions may also be indistinguishable, which would require further (or 
improved) measurements. Therefore the conditions needed to remove the exact degeneracies 
are necessary but not sufficient for removing degeneracies in an actual experiment. Further- 
more, parameter degeneracies occur nearly everywhere in the (5, #13) plane, since in most 
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TABLE I. Four classes of degenerate solutions which could potentially remain after three ap- 



pearance measurements. 



Case 


sgn(<5m| 1 ) 


02'J, 


I 


+ 


< 7T/ '4 


11 




< vr/ 4 


ill 


+ 


> 7T/ '4 


IV 




> tt/4 



cases any given ellipse (for fixed 13 ) overlaps an adjacent ellipse (with different 9 13) 



B. Parameter degeneracies with three measurements 

Since two measurements will generally have parameter degeneracies, the natural question 
is whether a third measurement sufficient to remove all of them. To extend the degeneracy 
analysis to more than two appearance measurements, what is needed is a straightforward 
and systematic method for determining whether or not a given set of such measurements 
resolves all of the degeneracies over the entire region of interest in the parameters S and #13. 
The method we present here uses the approximate expressions for the probabilities given in 
Sec. 0. 

We begin by analyzing the case of three appearance measurements. From the approxi- 
mate probability expressions we have three linear equations involving x 2 , x cos 5, and x sin 5. 
As long as the three equations are linearly independent (which can be realized if no two 
measurements simultaneously have the same L, E u , and type of neutrino) these are easily 
solved, giving a unique solution for x and 5. However, in principle four possible solutions 
may still exist, corresponding to the sgn(5m| 1 ) and #23 ambiguities (see Table D; eacn such 
solution gives identical predictions for the three probabilities being measured, but can have 
different values for 9 is and/or 5. 

To determine if the remaining degeneracies have been resolved by the third measurement, 
we use the fact that the variables x 2 , xcos5, and xsin.5 are not independent, and in the 
"true" solution they must obey the relation 

x 2 = (x cos 5) 2 + (x sin 5) 2 . (16) 

The question is then whether the "fake" solutions can also satisfy the constraint in Eq. (|16|). 



We have found that for many values of 5 and #13, the fake solutions do not obey Eq. (16) 
(thus removing the degeneracy). However, for some particular values of 5 and # 13 there 
exist fake solutions that also obey Eq. flTBD, so that degeneracy may remain even after three 
measurements; the values of 5 and #13 for which degenerate solutions still exist form lines 
in the (S, 6^3) plane. It is easy to understand why this is so: with two measurements, most 
points in the (5, 6*13) plane have degeneracies, but the condition in Eq. ([if]) imposes an 
additional constraint so that for three measurements the dimension of the degenerate space 
is reduced from two (the plane) to one (lines). Since the forms of P and P are the same 
(i.e., a linear combination of x 2 , xcos5, and xsin5), any combination of three neutrino 
and/or antineutrino appearance measurements will yield lines in (5, #13) space where the 
fake degenerate solutions are not eliminated by three measurements. 
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For example, consider neutrino measurements made at L = 730 km and E v = 2.66, 1.77, 
and 1.34 GeV (A = n/3, ir/2, and 27r/3, respectively); Fig. ||a shows the lines in (5, 9 13 ) 
space where a fake solution remains degenerate with the true solution (taken to be Case I 
in the figure). Figure shows similar curves for L = 730 km and E u = 3.54, 1.77, and 
0.89 GeV (A = 7r/4, 7r/2, and tt, respectively). In Fig. we assume 



\5m z 31 \ = 3 x 10~ 6 eV z , sin^ 29 23 = 0.90 , 
Sm 2 21 = 5 x 10" 5 eV 2 , sin 2 29 12 = 0.80 , (17) 

which allows the possibility of a 623 ambiguity as well as the (5, (9 13 ) and sgn(5m 2 1 ) ambi- 
guities; this will be the standard parameter set used in this paper unless stated otherwise. 
The key to determining if three measurements are sufficient to remove degeneracies in a 
particular experiment is whether or not the experiment is sensitive to the values of S and 
#13 where the lines of degeneracies occur. We see that in the two examples in Fig. |2|, degen- 
eracies remain for a wide range of 5 when sin 2 2# 13 > 0.01. Thus either of these two sets of 
measurements are not sufficient to ensure there are no degeneracies in the next generation 
of long-baseline neutrino experiments, which should be sensitive down to sin 2 29\3 = 0.01 
or below praMB0|. 



C. Parameter degeneracies with four measurements 

When three measurements are not sufficient to completely resolve the degeneracies, a 
fourth measurement must be made. A straightforward way to determine whether or not the 
fourth measurement has removed the remaining degeneracies is as follows. Take all possi- 
ble combinations of three measurements within the four (there are four such combinations) 
and perform the above analysis on each combination of three measurements. Each three- 
measurement subgroup may have values of 5 and 813 for which degeneracies remain (if not, 
then that subgroup alone is sufficient to remove the degeneracies, and a fourth measure- 
ment is not needed). If the degenerate lines in (<5, $13) space for the four subgroups do not 
intersect at a common point, then in principle the four measurements have resolved all the 
degeneracies; if not, then some degeneracies remain. 

For example, in the three measurements used to obtain Fig. |2]a there were degeneracies 
for a range of 5 when sin 2 2#i3 > 0.01. After a fourth measurement with L = 730 km 
and E v = 5.31 GeV (A = ir/Q) is added to the measurements used in Fig. ^|a, the lines of 
degeneracy for the four subgroups of three measurements are shown in Fig. |3|. There is a 
point near sin 2 2# 13 = 0.0085 and S = 0.197T where the lines of degeneracy all intersect when 
comparing I and IV, so that these two cases remain degenerate after these four measurements 
(see Fig. |||c), Furthermore, there is also a point near sin 2 2#i3 = 0.0015 and 5 = 0.497T where 
the four lines nearly intersect at the same point when comparing I and II, so that these two 
cases would be degenerate once experimental uncertainties are considered. 

It is not hard to understand why some degenerate points remain in (5,6*13) space after 
four measurements. The fourth data point adds one additional constraint on the parameters, 
which then reduces the dimension of the space that is degenerate from one (lines) to zero 
(points). For an experiment to be assured of resolving the parameter degeneracies, these 
degenerate points must lie at sin 2 29 13 below the experimentally accessible region. 
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We note that to determine the parameter degeneracies between, say, Case I and Case II 
with more than two measurements using a purely numerical approach, for each point in 
(5, 6*i 3 ) space for Case I, one would have to search the entire (5, 6*i 3 ) space for Case II 
to check for an identical set of probabilities. With our analytic approach, one need only 
search the parameter space once, since the existence of parameter degeneracies is determined 
algebraically. 

Figure [|a shows some typical parameter degeneracies between Cases I and II if only 
one measurement is made for v and v at L = 300 km and A = tx/2 (corresponding to 
E v = 0.73 GeV), assuming the same oscillation parameters as Eq. (|T7]) except #23 = 7r/4. 
We note that these degeneracies can mix CPV and CPC solutions, especially when 8 13 is not 
small; for some sets of parameters, CPC can give identical results to maximal CPV with 
the opposite sgn(5m| 1 ). To illustrate how additional measurements reduce the possibilities 
for degeneracies, Fig. f|b shows the parameter regions in Case I that have degeneracy with 
Case II, for two, three, and four measurements; the degenerate space is reduced from an area 
to lines to points. The corresponding degenerate parameters for Case II in Fig. |]b occur at 
the same #13 with 5 —>■ 5 ± 7r; it is not hard to show that this symmetrical situation is an 
artifact of choosing the same energy for v and and would not necessarily be true if the v 
and v energies were different. 

In summary, two measurements necessarily have parameter degeneracies over a broad 
range of (5, 9 13 ) space; adding a third measurement reduces the degeneracy to lines in the 
parameter space, and adding a fourth measurement reduces it further to isolated points. 
This behavior is to be expected since each additional measurement adds one more constraint, 
which has the effect of reducing the dimension of the degenerate subspace by one. Adding a 
fifth measurement with a linearly independent combination of x 2 , xcosS and xsinfi should 
then remove the degeneracy completely for all values of S and #13. For example, in Fig. |3|, if a 
fifth measurement at L = 730 km is made at, say, A = 5tt/6, then no point in the parameter 
space remains degenerate. Of course, practically speaking, less than five measurements may 
be sufficient to eliminate the possibility of parameter degeneracies in any given experiment if 
the degeneracies occur in a region of parameter space in which the experiment is not sensitive. 
In the next section we explore some specific cases with less than five measurements. 

IV. DETECTOR SCENARIOS FOR RESOLVING PARAMETER 

DEGENERACIES 

In the previous section we presented a general method for determining if a set of more 
than two measurements are sufficient to remove all of the parameter degeneracies. In this 
section we discuss specific detector scenarios using this method. We emphasize scenarios 
where at least some of the measurements can be made at the same L but different E u , since 
these are well-suited to off-axis beams (where different energies can be obtained by having 
a detector at different locations with respect to the beam axis). If a detector cluster could 
be used, the time required to make all of the necessary measurements could be reduced. We 
also compare the detector cluster scenarios to scenarios in which a detector cluster is not 
used. 
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A. Two v and two v measurements at one L 

The standard means for the study of CP violation is to measure both —>■ v e and 
Vn — > v e at the same L and E u . As we have seen, there could be as much as an eight-fold 
degeneracy in the parameters with just these two measurements, and it may be impossible 
to establish CPV from v and v measurements at a single L and E v . The simplest way one 
could imagine obtaining another set of independent measurements is to measure — > v e 
and/or — > v e at the same L but with different E u ; the detector and beamline remain 
the same, and only the beam energy and/or primary particle in the beam (z/ M or v^) would 
need to be changed from one run to the next. With a single detector this would take four 
separate runs. However, with two detectors at different locations with respect to the beam 
axis (so that each receives neutrinos of different energy) these four measurements could be 
taken with only two runs (one with v and one with v). 

Figure |5| shows the points in (5, (9 13 ) space where degeneracies remain due to the sgn^m^) 
and 823 ambiguities after measurements at two energies for both neutrinos and antineutrinos, 
assuming that v energies are the same as the v energies and that the measurements are all 
made at approximately the same L, for different choices of L. We examined scenarios 
where one energy was chosen so that Ai = tt/2, and the second energy was chosen so that 
A 2 = 7r/6, 7r/4, 7r/3, 27r/3, or it; some representative examples are shown in the figure. 
We chose Ai = tt/2 because the probabilities are large there (see / and / in Fig. |l|), the 
cos S term in the probabilities vanishes (thereby isolating the sin 5 term), and the survival 
channel —>■ is near an oscillation minimum. All possible degeneracies are included in 
Fig. H (e.g., I vs. II, II vs. I, I vs. Ill, etc., where the first case is the true solution), and 
different types of degeneracies are not distinguished. In all cases there are some points in 
the parameter space where degeneracies are still possible. We found that if A 2 < 7r/2, the 
degeneracies tended to be at smaller values of #13 than for A 2 > tt/2; the probabilities are 
also higher there (see / and / in Fig. [I]). 

For off-axis beams, both the neutrino flux and energy fall off with increasing off-axis 
angle (see Eqs. 1 and 2 of Ref. p3|); after eliminating the angular dependence one obtains 
for the energy dependence of the flux 

$oc££, (18) 

so that $ oc A~ 2 for measurements at the same L, Therefore if two off-axis measurements 
had very different A values, one would have a much reduced flux compared to the other. 
Hence it may be preferable if the ratio between the two values of A is not too large, since 
then the fluxes would be more equal. For example, with Ai — tt/2, A 2 = tt/4 and 7r/3 give 
very similar parameter degeneracies, but choosing A 2 = it/ 3 would make the fluxes of the 
two measurements closer; even then, the ratio of the fluxes would be 4:9. Of course, the 
detector further off axis could be made larger to help equalize the event rate. The detectors 
should not be too close in off-axis angle, either, since then the A values would be very 
similar and the expressions for the probabilities would no longer be linearly independent 
after experimental uncertainties are taken into account. 

Similar to the discussion of Fig. [|b, since the v energies are the same as the v energies, 
the degenerate points in Fig. [| come in symmetrical pairs, each with the same #13 and S 
differing by 7r. It can be shown that the degenerate parameters for Case I (III) when it is 
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degenerate with Case II (IV) exhibit this symmetry property with the parameters for II (IV) 
when it is degenerate with I (III). Therefore for degeneracies between I and II, or between 
III and IV, CPC and CPV solutions are not mixed, since 5 and 5 ± tt are either both 
CPC or CPV (although of course sgn^m^) is not determined). On the other hand, the 
symmetry also exists between the parameters for I when it is degenerate with IV (III) and 
the parameters for II when it is degenerate with III (IV), and between the parameters for IV 
(III) when it is degenerate with I and the parameters for III (IV) when it is degenerate with 
II; in all of these situations, the degenerate points that exhibit the symmetry property do not 
occur between two cases that are degenerate with each other, so there is still a CPC /CPV 
confusion for these degeneracies. 

The values of 6 at which degeneracies occur in Fig. || are insensitive to the assumed 
values of Sm^ and sin 2 2#i 2 . However, we found that the values of #13 at which degeneracies 
occur vary approximately as 



26fr n = (sin 2 20&°») o - T 2 , 1 T72 > (19) 




sm' 

where the degenerate parameter is (sin 2 28f^ 9en ) when Sm^i = 5 x 10~ 5 eV 2 and sin 2 2# 12 = 
0.80. Thus the values of 6*13 at which degeneracies occur is somewhat sensitive to sin 2 2# 12 
(which varies from 0.67 to 0.92 at the 2cr level ||) and very sensitive to Srn^i (which varies 
from 3 x 10 -5 to 2 x 10~ 4 eV 2 at the 2a level ||). We also found that a variation in dm^ 
of 10% (the expected precision of MINOS, ICARUS, and OPERA) caused sin 2 2#i 3 of the 
degenerate points to vary by 30% or less; smaller 5m 2 x moved the degeneracies to higher 
sin 2 2# 13 ; the largest variations occurred at lower sin 2 2# 13 . The corresponding variation in 
S was much less, of order a few percent. 

In Fig. H], all of the examples at shorter L have degeneracies for sin 2 29i3 ~ 0.01 or larger, 
which is within the expected range of future long-baseline experiments with superbeams. 
However, at L — 300 km with A = n, the degeneracy for larger # 13 occurs for sin 2 26* 13 > 0.1 
(around 0.13), above the current upper bound. The only degeneracies with sin 2 2#i3 < 0.1 in 
the L = 300 km case occur at much smaller 6*13 (sin 2 2# 13 ~ 0.0013; see Fig. |]c), just below 
the region accessible at the 3cr level in the Super JHF to Hyper-K experiment [p^] . Thus 
v and v measurements with both A = 7r/2 and tt at L = 300 km do not have degeneracy 
problems over the entire search region. However, the values of |/| and |/| are significantly 
smaller for A = 71 (by more than a factor of 10 compared to A = 7r/2; see Fig. [I]), so that the 
probabilities are also significantly smaller there, leading to significantly reduced statistics 
that will limit the sin 2 2#i3 reach of such an experiment. Furthermore, for slightly smaller 
values of <5m 21 or sin 2 26*i 2 the degeneracy at sin 2 29i 3 = 0.13 will be pushed below 0.10, into 
the search region. 

Since the matter effect increases with increasing L and 6 I 13 , one would expect that at 
larger L degenerate solutions with opposite sgn(5mg 1 ) can be differentiated for smaller values 
of #13 . This property is evident in Fig. Q e.g., for L = 2900 km degeneracies occur only 
at sin 2 2#i3 < 0.002-0.003. If such an experiment were sensitive only to sin 2 2#i3 above this 
value, it would not have a degeneracy problem. The sin 2 26*13 sensitivity depends on the 
total flux (which varies with L and off-axis angle), so a detailed calculation must be made 
to determine the existence of degeneracies in any given experimental situation. 
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There may also be approximate degeneracies that cannot be differentiated due to statisti- 
cal and systematic uncertainties. Here, too, one also needs to include details of a particular 
experimental situation to determine if approximate degeneracies are a problem; such an 
analysis is beyond the scope of this paper. However, to give an indication of the effect of 
approximate degeneracies, in Fig. |6| we show regions in (5, #13) space where the differences 
between (x cos 5) 2 + (x sin S) 2 and x 2 for the false solution, when added in quadrature, are less 
than 10%. We see that in Figs. |^a and ^|b the approximate regions extend to sin 2 2813 that 
are nearly twice as large as the maximum values for the exact degeneracies. Furthermore, in 
Fig. |5|c there are no exact degeneracies in the region shown, but approximate degeneracies 
extend to sin 2 28 13 as large as 0.007. Therefore it is clear that approximate degeneracies will 
occur at somewhat higher #13 than the exact degeneracies. 



B. Four v measurements at one L 

Antineutrino measurements have two drawbacks compared to neutrino measurements: 
the v cross section is approximately half as large, and v beam fluxes are less. Thus to achieve 
comparable statistics the running time for an antineutrino measurement must be more than 
twice as long as a similar neutrino measurement (apart from the possible enhancement of 
antineutrino oscillation probabilities at long L if Sm^ < 0). If it were possible to remove 
the degeneracies with neutrino measurements alone, the total time to complete all of the 
necessary measurements would be substantially reduced. With a single detector, this would 
require more than two separate runs (since, as we have shown, any two measurements alone 
will have degeneracies), each at a different E v . With two detectors at different locations with 
respect to the beam axis, only two runs would give four independent measurements, and only 
the on-axis beam energy or the detector off-axis angles would have to be changed between 
runs (although the latter might prove to be unfeasible). Another interesting possibility is to 
have a four-detector cluster, so that only a single neutrino run would be needed to eliminate 
all degeneracies. Although this many detectors might seem to be economically unfeasible, 
the running time is shorter and the beam would not have to be reconfigured, which would 
offset increased detector costs. In this section we examine the ability of four v runs to resolve 
degeneracies. 

In Fig. [7| we show the points in the (S, 813) plane that have degeneracy with four neutrino 
measurements at the same L, for several choices of L. We see that in all cases there are 
points in the parameter space that have degeneracies even after four measurements are 
taken. Similar to the case with two v and two v measurements, the degenerate solutions are 
pushed to lower 813 as L is increased, due to the larger matter effect, which helps distinguish 
sgn(5m| 1 ). For L = 2900 km, in the best cases we found that the largest sin 2 2# 13 with 
a degeneracy is 0.004-0.006. A comparison of Figs. ^] and shows that two v and two v 
measurements at the same L generally has degeneracies occurring at lower sin 2 2#i3 than 
four v measurements at one L. 

For a small off-axis angle 8, the neutrino energy is given approximately by |23 



E » - i 1 , Pfl /L , 2 , (20) 



0.43^ 

where E v is the pion energy in the lab frame and m n is the pion mass. Four different 
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measurements can be made simultaneously if four detectors are placed at different off-axis 
angles, and it is always possible to choose E n and the angles to select any desired combination 
of Aj. This is not the case if only two detectors are used and two runs are made with different 
E n (assuming the detectors stay at the same off-axis angles in both runs). From Eq. (f20[), 
the condition that a set of four Aj can be obtained in a two-detector, two-run scenario is 

A 1 (A 4 -A 1 )> A 2 (A 3 -A 2 ), (21) 

where A 4 > A 3 > A 2 > A 1 . In practice, Eq. (^) is not hard to satisfy, especially if A 4 is 
not too small and the difference between A 3 and A 2 is not too large. All of our examples in 
Fig. [7| can be obtained in two runs with two detectors. 

One drawback of making only neutrino measurements is that this is not a direct mea- 
surement of CPV (such as would be the case of having both v and v measurements). 
Furthermore, if Sm^ < 0, then all of the measurements would have a suppression due to 
matter effects, especially at longer L. A detailed calculation including experimental uncer- 
tainties would have to be made to see if the improved resolution of degeneracies at longer L 
can overcome the rate loss if Sm^ < 0. 

C. Measurements at more than one L 

In the previous examples all of the detectors were assumed to be at the same distance. 
As we have seen, measurements at longer L push the degeneracies to lower sin 2 2#i 3 . Here we 
consider having two v measurements made at one distance (simultaneously with a detector 
cluster) and a third v measurement at another distance. We examine two possibilities: one 
with the detector cluster at the shorter distance, the other with the detector cluster at the 
longer distance. 

In all cases we assume the first measurement is made with a v beam at L\ = 300 km and 
Ai = 7r/2 (the approximate values for the proposed JHF to Super-K experiment). If the 
second v measurement is also at L 2 = 300 km with A 2 ^ A 1; and a third v measurement at 
L 3 > L 2 , we found that the degeneracies are less severe (i.e., they occur at smaller sin 2 26> 13 ) 
when A 3 = A 2 < 7r/2; some of the better examples we found are shown in Table f|. For the 
longer distance we used the L values considered in Fig. [7|, plus 2100 km (nominally JHF to 
Beijing). Having A 2 < n/2 and A 3 = n/2 also gave good results. In the cases we tested 
with the second and third v measurements both at longer L, the best results were achieved 
if one of the measurements at the longer distance had A = 7r/2; some of the better examples 
we found are shown in Table [TTT]. 

In both Tables |T| and |T|, experimental L and E v choices can be made such that degen- 
eracies occur only at sin 2 26*i 3 = 0.003 or less (assuming standard values for sin 2 2#i 2 and 
5m 2 2 ). Thus, although there are only three measurements, with appropriate choices of L and 
A, degeneracies are less of a problem (i.e., they occur at smaller # 13 ) than for the scenarios 
where all measurements are done at the same L. Of course, for only three measurements the 
regions where degeneracies occur are lines in (5, #i 3 ) space, rather than isolated points, so 
they are more likely to occur in experiments with sin 2 2#i 3 sensitivity of order 0.003 or less. 
Furthermore, as with the scenario with four v measurements at one L, a determination of 
CPV is indirect (depending on the three-neutrino parametrization) and probabilities would 
be suppressed if < 0. 
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TABLE II. Maximum sin 2 2#i3 that has a degeneracy when the first v measurement is made 
at L\ = 300 km and Ai = vr/2, a second at L2 = 300 km with A2 7^ Ai , and a third at L3 > L2, 
for different choices of the second and third measurements. The other neutrino parameters are the 



same as in Eq. (p~7|) . 











L 3 








AsM 


730 km 


1290 km 


1770 km 


2100 km 


2900 km 


tt/4 


tt/4 


.0012 


.0013 


.0013 


.0013 


.0013 




tt/2 


.0028 


.0033 


.0038 


.0042 


.0060 


tt/3 


tt/3 


.0011 


.0012 


.0012 


.0013 


.0013 




tt/2 


.0021 


.0025 


.0030 


.0034 


.0047 



TABLE III. Maximum sin 2 2#i3 that has a degeneracy when the first v measurement is made 
at L\ = 300 km and Ai = vr/2, and two other v measurements are made at L2 = L3 > L\, for 
different choices of A for the second and third measurements. The other neutrino parameters are 
the same as in Eq. (|i~7|). 



A2H 


A 3 H 


730 km 


1290 km 


L2 = Lz 
1770 km 


2100 km 


2900 km 


tt/2 


tt/4 


.0027 


.0030 


.0033 


.0038 


.0047 


tt/2 


tt/3 


.0021 


.0022 


.0025 


.0027 


.0035 


tt/2 


2tt/3 


.0030 


.0042 


.0067 


.0094 


.032 



D. Scenarios without a detector cluster 

All of our previous examples included some measurements at more than one energy, but 
at the same L, which is particularly suited to detector cluster experiments. Here we examine 
scenarios which do not use a detector cluster. 

First we consider a v and v measurement at one L and E v , and another v measurement 
at the same L and different E,,. This is similar to the case of two v and two v measurements 



all at the same L discussed in Sec. [IV A| , except that the second v measurement is missing, 
and the two v measurements are done in separate runs with different v energies. Since 
there are only three measurements, the parameter degeneracies are lines in (S, #i 3 ) space. 
We examined cases where v and v measurements were done at A = tt/2 and another v 
measurement was done at A = rr/6, tt/4, tt/3, 2tt/3, 3tt/4, or rr, for L = 300, 730, 1290, 
1770, 2100, or 2900 km. We found that for L < 1290 km, in the best cases degeneracies 
can occur for sin 2 2# 13 as high as 0.01 (0.10 for L = 300 km). For the best cases with 
L > 1770 km degeneracies can occur for sin 2 2# 13 as high as 0.005-0.008, which is a factor 
of two or more higher than for two v and two v energies at the same L. Thus making a 
second v measurement not only reduces the set of degenerate points from lines to points, 
but it also pushes the degeneracies to lower sin 2 26 l 1 3. 

Next we consider a v and v measurement at one L and E U1 and another v measurement 
at a different L. The neutrino energies for the two v measurements may or may not be 
the same. We examine two possibilities, where the v measurement is done at either the 
shorter or longer L. In all cases we assume the first measurement is made with a v beam at 
L\ = 300 km and Ai = tt/2. If the v measurement is also at L2 = 300 km with A2 = 7r/2, 
and a v measurement at L 3 > 300 km, we found that for A 3 ~ tt/3 (A3 ~ 2tt/3) the 
degeneracies are less severe at larger (intermediate) L 3 ; see Table |V[ Another possibility 



is to have the v measurement done at the longer L. Typical results are shown in Table 0. 
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TABLE IV. Maximum sin 2 2#i3 that has a degeneracy when v and v measurements are made 
at L\ = L2 = 300 km and Ai = A2 = vr/2 and another v measurement is made at L3 > 300 km, 
for different choices of L3 and A3. The other neutrino parameters are the same as in Eq. (|17|). 









Lz 








730 km 


1290 km 


1770 km 


2100 km 


2900 km 


tt/3 


.089 


.018 


.008 


.006 


.004 


2vr/3 


.035 


.010 


.016 


.025 


.071 



TABLE V. Maximum sin 2 29\3 that has a degeneracy when a v measurement is made at 
L\ = 300 km and Ai = n/2, and v and v measurements are made at L2 = L3 > 300 km, for 
different choices of the second and third measurements. The other neutrino parameters are the 
same as in Eq. ([l7|). 











L2 = Lz 








A3 (p) 


730 km 


1290 km 


1770 km 


2100 km 


2900 km 


7T/3 


tt/2 


.050 


.014 


.005 


.006 


.010 


2tt/3 


tt/2 


.038 


.010 


.016 


.022 


.067 


tt/4 


tt/4 


.056 


.013 


.007 


.005 


.008 


tt/3 


tt/3 


.038 


.008 


.005 


.006 


.013 


2tt/3 


2tt/3 


.033 


.011 


.011 


.015 


.063 



Generally speaking the measurements in Tables [IV] and [V| have degeneracies for higher 



sin 28 13 than with three v measurements, one at a different L (Tables [TT] and |TT1|). 



V. SUMMARY AND DISCUSSION 



We summarize the important points of our paper as follows: 

(i) For any two appearance measurements there may be as much as an eight-fold param- 
eter degeneracy (resulting from simultaneous (5, #13), sgn(5m 2 1 ) and (6*23, 7r/2 — # 23 ) 
ambiguities) for most points in the (5, #13) plane. Making a third appearance mea- 
surement resolves the (5, #13) ambiguity and reduces the regions where the remaining 
degeneracies occur to lines in (5, #13) space. Making a fourth appearance measurement 
reduces these degeneracies to isolated points in the (8, #13) plane; a fifth measurement 
then in principle removes all remaining degeneracies. 

(ii) Two v and two v measurements at the same L can be made so that parameter degenera- 
cies only occur at isolated points in (5, 813) space at sin 2 26>i 3 < 0.01-0.02 (0.002-0.003) 
for L = 300 (2900) km. If the v energies are the same as the v energies, there will be 
no CPC/CPV confusion for degeneracies between Cases I and II (denned in Table |), 
or between Cases III and IV, although sgn(5m 2 1 ) is not determined. This scenario 
could be completed in two runs (one with a v beam and one with a v beam) with two 
detectors at different positions with respect to the beam axis (a two-detector cluster). 

(iii) Four v measurements at the same L can be made so that parameter degeneracies 
only occur at isolated points in (5, # 13 ) space at sin 2 2# 13 < 0.04 (0.004) for L = 
300 (2900) km. Such a scenario could be implemented in two runs (two v beams with 
different on-axis energies) with a two-detector cluster, or in one run (a single v beam) 
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with a four-detector cluster. It would have the advantage that the v cross section and 
flux would be larger than for v\ disadvantages include the fact that it is an indirect 
measurement of CPV and it would have rates suppressed by matter effects at longer 
L if <5m^ < 0. 

(iv) Two v measurements at one L and a third at another L can be made so that pa- 
rameter degeneracies only occur at sin 2 2^13 < 0.001-0.003. Since there are three 
measurements, the degeneracies occur along lines in the (S, 813) plane. Such a scenario 
could be implemented in two runs, with a two-detector cluster at one L and a single 
detector at the other L. This scenario uses the matter effect at longer L to help push 
degeneracies to lower #13 , and, since it uses only v beams, has similar advantages and 
disadvantages as (iii). 

(v) Both a v oscillation measurement and a v oscillation measurement at one L and a v 
measurement at another L can be made so that parameter degeneracies occur along 
lines in the (5, #13) plane for sin 2 2# 13 < 0.005-0.010. Similar results can be obtained for 
two separate v measurements and one v measurement at the same L. These scenarios 
do not use a detector cluster. 

(vi) All of the examples we show assume 8rrv^ x = 5 x 10~ 5 eV 2 and sin 2 2# 12 = 0.80. We 
found that the maximum # 13 that may have parameter degeneracies varied strongly 
with Sm^, and less so with sin 2 2# 12 (see Eq. [L9]); the corresponding 5 values were 
unaffected. There was also a dependence on Sm 3l , but if is known to 10% its 
uncertainty does not greatly affect our degeneracy analysis. 



A summary of the requirements for the scenarios discussed in Sec. |TV|is given in Table [VI|; 
also shown is the best-case (i.e., lowest) maximum value of sin 2 2#i3 that we found for which 
a parameter degeneracy could occur, assuming the standard parameter set of Eq. (|17|). 
Generally speaking, we found that measurements at more than one L did better at resolving 
parameter degeneracies than if all measurements were done at a single L; for example, a 
third v measurement at a different L than the first two measurements appears to be better 
than a third and a fourth at the same L. Also, scenarios with detector clusters resolved 
degeneracies as well as or better than those without a detector cluster, and required fewer 
runs. Most scenarios resolved parameter degeneracies the best at longer L (1770 km or 
more), except for the scenario with a two-detector cluster at one L and a single detector at 
another L, which actually did better when the second distance was not as long (1290 km or 
less) . 

If 6*23 — 7r/4 (the value favored by atmospheric neutrino experiments), then the 823 
ambiguity vanishes and the only two cases to consider after three measurements are I and II 
(the sgn(5m 2 1 ) ambiguity). We have also performed the degeneracy analyses for sin 2 2# 23 = 
1, and found that all of the relevant degeneracies are at least a factor of two lower in sin 2 28 13 



than shown in Table |V| for the shorter distances (L < 1290 km), and in most cases lower 
than that for the longer distances. Therefore if it is known that 823 is close to it/ 4 (for 
example, from measurements at MINOS or ICARUS), the parameter degeneracy problem is 
greatly reduced. 

Some of the scenarios we discussed involve measurements with neutrinos only (not an- 
tineutrinos). Although it is possible to extract the neutrino mixing angles and CP phase 
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TABLE VI. A summary of detector scenarios discussed in this paper. The right-most col- 
umn indicates the best-case (i.e., lowest) maximum value of sin 2 2^13 that we found for which a 



parameter degeneracy could occur, assuming the parameters of Eq. fllTT) . 













■ 2 n/l 

max. sin zt/13 


# of 


# of 


# of 


# of 


# of 


without 


detectors 


beams 


beamlines 


runs 


measurements 


degeneracy 


2 at one L 


2 (v and v) 


1 


2 


4 


.002 


2 at one L 


2 (1/ at Ei and £? 2 ) 


1 


2 


4 


.004 


4 at one L 


1(1/) 


1 


1 


4 


.004 


2 at Li, 1 at L2 


1(1/) 


2 


2 


3 


.001 


1 at one L 


3 (1/ at £1 and £"2, i>) 


1 


3 


3 


.005 


1 at Li, 1 at L2 


2 (z/ and p) or 3 (2 u and 1 i>) 


2 


3 


3 


.004 



from such measurements, the result relies on the three-neutrino mixing assumption. A 
more robust determination of CP violation which would not be as model dependent would 
include a measurement involving antineutrinos; then the measurement of an asymmetry 
between neutrino and antineutrino oscillation probabilities (after correcting for the CP vi- 
olation induced by matter effects) would be direct and definitive. 

The detector cluster scenarios are especially well-suited for detector designs that em- 
phasize large, cheaply-built detectors. We note that the detectors in a cluster would not 
necessarily have to be the same size; detectors at larger off-axis angles could be made larger 
to compensate for the reduced flux off-axis. 

We emphasize that our analysis only considered exact degeneracies for monoenergetic 
beams. Experimental uncertainties will expand the regions with degeneracies, especially 
at small sin 2 26^ (where probabilities and hence events rates are lower) and at longer L 
(where beam fluxes fall off). On the other hand, neutrino beams are not precisely mo- 
noenergetic, and energy spectrum information can help resolve the degeneracies, even for 
relatively narrow spectrum beams [ID]; the survival channel — > can in principle also 



provide additional information, although the effects of #13 and S are not at leading order. 
A definitive analysis would have to include all of these factors to determine if parameter 
degeneracies are a problem in a particular experiment. The results of this paper can serve 
as a guideline for which detector/beam scenarios are likely to encounter difficulties with 
parameter degeneracies. 



ACKNOWLEDGMENTS 



This research was supported in part by the U.S. Department of Energy under Grants 
No. DE-FG02-95ER40896, No. DE-FG02-01ER41155 and No. DE-FG02-91ER40676, and 
in part by the University of Wisconsin Research Committee with funds granted by the 
Wisconsin Alumni Research Foundation. 



17 



REFERENCES 



[1] Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 87, 071301 (2001) ||arXiv:nucl- 
ex/0l060l5|. 



[2] Q. R. Ahmad et al. [SNO Collaboration], |arXiv:nucl-ex/0204008| , |arXiv:nucl- 
ex/0204009l 



[3] J. N. Bahcall, M. H. Pinsonneault and S. Basu, Astrophys. J. 555, 990 (2001) 

||arXiv:astro-ph/0010346| . 
[4] V. Barger, D. Marfatia, K. Whisnant, and B.P. Wood, |arXiv: hep-ph / 0204253| . 
[5] A. Bandyopadhyay, S. Choubey, S. Goswami and D.P. Roy, |arXiv:hep-ph/0204286| ; 

J. N. Bahcall, M. C. Gonzalez-Garcia and C. Pena-Garay, arXiv:hcp-ph/0204314 ; 

P. Aliani, V. Antonelli, R. Ferrari, M. Picariello, and E. Torrente-Lujan, |arXiv:he~j> 

ph/0205053| ; PC. de Holanda and A. Yu. Smirnov, |arXiv:hep-ph/0205241| . 
[6] P. Alivisatos et al, STANFORD-HEP-98-03; V. Barger, D. Marfatia and B. P. Wood, 

Phys. Lett. B 498, 53 (2001) | |arXiv:hep-ph/00ll251| . 
[7] T. Toshito [Super-Kamiokande Collaboration] , |arXiv: hep-ex /0105023 . 



J. E. Hill [K2K Collaboration], |arXiv: hep-ex/ 110034 



[9] M. Apollonio et al. [CHOOZ Collaboration], Phys. Lett. B 466, 415 (1999) ||arXiv:hep- 



ex/9907037 



[10] F. Boehm et al, Phys. Rev. D 64, 112001 (2001) ||arXiv:hep-ex/0107009i . 
[11] MINOS Collaboration, Fermilab Report No. NuMI-L-375 (1998). 

[12] A. Rubbia for the ICARUS Collaboration, talk given at Skandinavian NeutrinO Work- 
shop (SNOW), Uppsala, Sweden, February (2001), which are available at [http:/ /pcno- 
| meth4.cern.ch/publicpdf.html] . 

[13] OPERA Collaboration, CERN/SPSC 2000-028, SPSC/P318, LNGS P25/2000, July, 
2000. 

[14] V. Barger, A. M. Gago, D. Marfatia, W. J. Teves, B. P. Wood and R. Z. Funchal, 
larXiv: hep-ph /0 1 1 03"93| . 



[15] J. Burguet-Castell, M. B. Gavela, J. J. Gomez-Cadenas, P. Hernandez and O. Mena, 
Nucl. Phys. B 608, 301 (2001) [|arXiv:hep-ph/0103258 . 



[16] V. Barger, D. Marfatia and K. Whisnant, |arXiv:hep-ph/010809U . 

[17] V. Barger, D. Marfatia, and K. Whisnant, Phys. Rev. D 65, 073023 (2002) ||arXiv:hep- 

| ph/0112111 . 

[18] T. Kajita, H. Minakata, and H. Nunokawa, |arXiv: hep-ph/ 112345 . 

[19] P. Huber, M. Lindner, and W. Winter, |arXiv:hep-ph / 0204354 

[20] H. Minakata and H. Nunokawa, JHEP 0110, 001 (2001) ||arXiv:hep-ph/0108085|| . 

[21] D. Beavis et al, E889 Collaboration, BNL preprint BNL-52459, April 1995. 

[22] Y. Itow et al, |arXiv:hep-ex/0106019| . 

[23] A. Para and M. Szleper, |arXiv:hep-ex/0110032 . 



[24] G. Barenboim, A. de Gouvea, M. Szleper, and M. Velasco, |arXiv:hep-ph/ 0204208 



[25] M. Aoki, K. Hagiwara, Y. Hayoto, T. Kobayashi, T. Nakaya, K. Nishikawa, and N. Oka- 

mura, |arXiv:hep-ph /01 12338| . 
[26] A. Cervera et al Nucl. Phys. B 579, 17 (2000) [Erratum-ibid. B 593, 731 (2000)] 

||arXiv:hep-ph/0002l08|| . 

[27] M. Freund, Phys. Rev. D 64, 053003 (2001) [|arXiv:hep-ph/010330q] . 

[28] A. Dziewonski and D. Anderson, Phys. Earth Planet. Inter. 25, 297 (1981). 



18 



[29] K. Kimura, A. Takamura, and H. Yokomakura, |arXiv:hep-ph/0205295 . 
[30] V. Barger, S. Geer, R. Raja, and K. Whisnant, Phys. Rev. D 63, 113011(2001) 
||arXiv:hep-ph/00120171 . 



19 




I vs. II - 
I vs. Ill 
I vs. IV 






sin 28 13 



sin 28,. 



FIG. 2. Lines in (sin 2 28 13 , 5) space for Case I where three neutrino measurements are not 
sufficient to eliminate the ambiguity between I and II (solid curves), I and III (dashed), and I and 
IV (dotted), assuming L = 730 km and that the other parameters have the values given in Eq. (|l7|). 
The values of E u are chosen so that the oscillation arguments are (a) A = tt/3, vr/2, and 2-7r/3, 
[E v = 2.66, 1.77 and 1.33 GeV, respectively), and (b) A = vr/4, vr/2, and ir (E u = 3.54, 1.77 and 
0.89 GeV, respectively). 
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FIG. 3. Lines in (sin 2 2#i3, 5) space for Case I where subgroups of three neutrino measurements 
are not sufficient to eliminate the parameter degeneracies between (a) I and II, (b) I and III, and 
(c) I and IV, assuming L = 730 km and that the other parameters are the same as in Eq. (|l7|). 
The values of E v for the four measurements are chosen so that the oscillation arguments are 
A = 7r/6, 7r/3, 7r/2, and 2-7r/3, (E u = 5.31, 2.66, 1.77 and 1.33 GeV, respectively), and the four 
subgroups are labeled as follows: E v = 5.31 GeV omitted (solid), E v = 2.66 GeV omitted (dashed), 
E v = 1.77 GeV omitted (dotted), E v = 1.33 GeV omitted (dash-dotted). 
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sin 2 26 13 sin 2 29 13 

FIG. 4. (a) Sample points in (<5, 6*13) space that have parameter degeneracies between Cases I 
(boxes) and II (crosses), defined in Table |, after v and v measurements at L = 300 km and A = it/2 
(corresponding to E u = 0.73 GeV), assuming the other parameters are the same as in Eq. (|i~7|), 
except for #23 = vr/4. The sets of parameters that are degenerate with each other are linked by a 
line. For each pair of linked points there is another pair with 5 — > ir — 5 (not shown), (b) Region 
in Case I parameter space that have degeneracies with Case II after (i) the two measurements in 
(a) (area between the solid curves), (ii) an additional v measurement at L = 300 km and A = tt/3 
(along dashed curves), and (iii) a final v measurement at L = 300 km and A = n/3 (boxes). The 
corresponding regions for Case II are found by taking 5 — > 8 ± ir. We note that after only the first 
two measurements there is also a (5,ir — 5) ambiguity everywhere in the plane, which is removed 
by the third measurement. 
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FIG. 5. Points in (sin 2 29 13 , 6) space where two v and two v measurements are not sufficient to 
eliminate all of the parameter degeneracies, for several values of L, assuming the other parameters 
are the same as in Eq. (|l7|). In each case measurements assumed to be made at two energies such 
that A = vr/2 and (a) A = vr/3, (b) A = 2vr/3, and (c) A = vr. 
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FIG. 6. Regions in (sin 2 2#i3, <5) space where two v and two v measurements at L = 300 km 
with A = 7r/2 and A = ir/3 still have approximate parameter degeneracies, assuming the other 
parameters are the same as in Eq. (|l7|). The three degeneracies shown are (a) I vs. II (I is the true 
solution and II is the false solution), (b) III vs. IV, and (c) III vs. II. The condition defining the 
approximate region is given in the text. Exact parameter degeneracies are indicated with pluses. 
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FIG. 7. Points in (sin 2 29^, 5) space where four v measurements are not sufficient to eliminate 
all of the parameter degeneracies, for several values of L, assuming the other parameters are the 
same as in Eq. (|i~7|). In each case measurements are made at four energies such that (a) A = 7r/4, 
7t/3, 7t/2, and 37r/4, (b) A = tt/6, n/3, 7r/2, and 2-7r/3, and (c) A = ir/G, n/3, tt/2, and 7r. 
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